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^| The quantum dynamics of a simplest dissipative system, a particle moving in 



a constant external field , is exactly studied by taking into account its interaction 



with a bath of Ohmic spectral density. We apply the main idea and methods de- 

veloped in our recent work [1] to quantum dissipative system with constant external 

field. Quantizing the dissipative system we obtain the simple and exact solutions for 

the coordinate operator of the system in Heisenberg picture and the wave function 

of the composite system of system and bath in Schroedinger picture. An effective 

Hamiltonian for the dissipative system is explicitly derived from these solutions with 

Heisenberg picture method and thereby the meaning of the wavefunction governed 
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by it is clarified by analyzing the effect of the Brownian motion. Especially, the 
general effective Hamiltonian for the case with arbitrary potential is directly derived 
with this method for the case when the Brownian motion can be ignored. Using this 
effective Hamiltonian, we show an interesting fact that the dissipation suppresses the 
wave packet spreading. 
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1. Introduction 



This paper is mainly devoted to the application and generalization of the idea 
and methods developed in our recent work [1] to a dissipative system with constant 
external field. Though most discussions in this paper are proceeded based on ref. [1] , 
the context is written in a self-contained form. In ref. [1], we work on the case of 
harmonic oscillator moving in a bath with the Ohmic spectral density presented by 
Caldeira-Leggett [2]. It has been shown there that the wavfunction of system plus 
bath is precisely described by a direct product of two independent Hilbert space, 
one of which is described by an effective Hamiltonian while the other represents 
the effect of the bath, i.e., the Brownian motion. Therefore, this study clarifies 
the structure of the wavefunction of the system dissipated by its interaction with 
the bath and thereby reveals the relationship between the different approaches for 
quantum dissipative systems before. No path integral technology is needed in this 
treatment. Notice that the study of dissipative quantum systems(DQS), especially 
for the damped harmonic quantum oscillator (DHQO), has a rather long history [3-18] 
and has been paid much attention more recently due to the work by Caldeira and 
Leggett [2]. 

Now let us briefly describe two main approaches for quantum dissipative systems 
before the work in ref.[l]. To reproduce and quantize the phenomenological dissipa- 
tive equation 

Mq(t) = -m(t)-^-, (1.1) 

for one-dimensional dissipative system S with coordinate q, mass M and potential 
V(q), one approach is to embed it into an environment, a bath B of N harmonic oscil- 
lators interacting with the system S through certain coupling [3-12]. Since the bath 
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and system constitute a closed composite system C equals to S plus B , the quanti- 
zation of C is naturally direct. Then, its corresponding Heisenberg equation results 
in the above phenomenological dissipative equation in the operator form by eliminat- 
ing the variables of the bath through proper approximations, such as the Markovian 
approximation and the Wigner-Wisskopf approximation [9]. This approach provides 
the motion equation with both the friction force in dissipation process and the fluc- 
tuation force in Brownian motion. The path integral technique [2,13] and field theory 
method [14] were used in this approach sometimes. Along this direction, the work of 
Caldeira and Leggett reveals a remarkable fact that the dissipation can occur exactly, 
instead of approximately, if the spectral density of the bath is Ohmic (to be described 
later). It enlightens people to reconsider many problems of dissipation in an exact 
version. 

Another approach, for DHQO, is the use of an effective Hamiltonian [15-17] 

H E = H E (t) = ^-^V + l -Mu 2 e^ M q 2 . (1.2) 

which is now called Caldirora-Kani (CK) Hamiltonian [13]. With the canonical com- 
mutation relation 

[q,p]=ih, (1.3) 

this Hamiltonian automatically yields the dissipation equation (1.1) through the 
Heisenberg equation. Notice that the alternative forms of the effective Hamiltonian 
have been by many authors [17] and an elegant example can be found in ref.[18]. 
Though this approach is very convenient to treat some dynamical problems of dissi- 
pation process for both classical and quantum cases, such as tunnelling and motion of 
wave packet , it is only a purely phenomenological method without the microscopic 
mechanism constructing H E and thus the Brownian motion can not be analysed by 
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making its use. Especially, the meaning of wavefunction of the evolution governed 
the CK Hamiltonian is much ambiguous. 

The present paper applies the method and idea developed in reference [1] to dis- 
sipative system with a constant external field, to establish the connection between 
the above two approaches and to consider how the interaction between the bath and 
the system leads to an explicit description for the dissipative system in terms of the 
effective Hamiltonian. All the discussions in this paper is proceeded with a simplest 
model, a particle moving in one dimension with constant potential field, but the main 
idea and methods can be generalized for other cases. The present work is completed 
in an unified framework in accompany with the answers to the entangling questions 
listed as follows: 

i. How to quantize the dissipative system with the effective Hamiltonian in an exact 
context? 

ii. What is the meaning of the wavefunction governed by the effective Hamiltonian ? 

iii. How to construct the propagator for the dissipative system in the Heisenberg 
representation ? 

iv . What happens to the spreading of the wave packet in present of dissipation ? 
2. Exact Classical Langevin Equation 

In this section we reformulate the first approach about quantum dissipative system 
in Ohmic case with a simplest example that the dissipative system S is considered as 
a charged particle with unit mass, unit negative charge and coordinate q in a constant 
electric field E. It interacts with a bath B of N harmonic oscillators Bi of coordinates 
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Xj, mass rrij and frequency ojj. Let p and pj be the corresponding momentums to q 
and Xj respectively. The Hamiltonian of the composite system of S plus B is written 

as 

1 N v 2 1 

H = -p 2 - Eq + g[J2L- + -m J cu 2 (x J - q) 2 ] 

1 N 
= -p 2 -Eq + H B -J2 C 3 x,q + AV. (2.1) 



where 



= V — ^g 2 = V -m,a; 2 g 2 



The problem of renormalization of interaction for N = oo mentioned in ref.[2] is 
naturally enjoyed by the frequency-dependent coupling constant Cj = rrijUJ^ and the 
renormalized potential 

N r*1 N -i 

in the Hamiltonian (2.1). 

Through the canonical equation, the Hamiltonian (2.1) defines a system of classical 
motion equations 

N N (J2 

S = £-E^'-E-7> (2.2a) 

j=l j=l m 3 U 3 

c 

Xj = -uj 2 Xj -q. (2.26) 

rrij 

Notice that the canonical equation for the closed composite system C of S plus B 
determines the ordinary momentum-velocity relations 

P(t)=q(t), Pj{t) =mjXj{t). 

By making use of the Laplace transformations of eqs.(2.2), a direct of substitution of 
eq.(2.2b) into eq.(2.2a) yields an exact motion equation for the system S 

q = E + L(q) + G(t), (2.3) 



where 

N x(0) 
G(t) = £ Cjlxj^cosujjt + -J^lsirujjt], (2.4) 



and the second term in the right hand side in eq.(2.3) 



N 



Cls 2 



^) = p-H 2 E m ^ + , |) g -Ml (2-5) 



is determined by the inverse of Laplace transformation p 

p[q(t)\ = q(s) = r q (t)e~ st dt 



Usually, for the finite N or the general spectral distribution p(uj) of infinite os- 
cillators in bath, the dissipative term — rjq for a positive number 77 does not appear 
exactly. However, according to Caldeira and Leggett, a specific spectral distribution 
of the bath 



which is called Ohmic distribution , enable the sum over index j 



1i mjUj](s 2 + u]) 
to become an integral 

r 1 \ ° 2 i s2 j 2r l s2 r duJ i 
Jo rrijUjjis 2 + uf) n Jo (ouj + s 2 ) 

and thus results in 

L(q) = -m - vq(o)S(t) 

immediately through an inverse Laplace transformation. Then, the dissipative equa- 
tion - the classical Langevin equation 

q = E- V q + G(t), (2.7) 
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occurs as an exact evolution based on the elementary Hamiltonian dynamics. Because 
the dissipative process is invertable, it is significant to pay our attention only to the 
process with t > as follows. In this sense the impact S(t) in the dissipative equation 
does not play role in dynamical problems. 

Notice that the above Ohmic distribution (2.6) is only an alternative, but explicit 
and convenient reformulation of the Caldeira-Leggett's constraint 

J{uj) = r\uj 

on spectral density 

N (J2 

J(u) = -'£-f-S(u-u j ). (2.8) 

It is also remarked that the fluctuating external force G(t) acting on system S is 
due to the effect of the bath and depends on the initial states of the oscillators in 
the bath. For the classical statistical thermal average < > classical, G(t) obeys the 
dissipation-fluctuation relation at temperature T 

< G(t)G(f) >dass i cai= \r]KT8{t - t') (2.9) 

and < G(t) > c iassicai= 0. It actually is the classical Brownian force in dissipation 
process. In concluding this section it is pointed out that the classical Langevin 
equation has been studied for a long time by many authors by certain approximations 
, but here our emphasis is that it can appear exactly with Ohmic spectral density for 
the bath. 

3. Quantization of Dissipation with Effective Hamiltonian 

In the following sections, we first detail the description in ref.[l] of the the system 
plus bath with the present example. Then, we develope this theory to derive a general 
effective Hamiltonian for arbitrary potential. 



It can be observed from eq.(2.7) that, in Ohmic case, the action of bath B on 
the system S can be exactly described as two parts, the dissipative term — r\q only 
depending on the state of S and the Brownian force G(t) depending on the initial 
state of bath. In order to study various dynamical problems in dissipative process, it 
is necessary to determine in what sense the dissipative system can be isolated from 
the environment as a 'quasi-closed' system only depending on its own variables. Such 
a 'quasi-closed' system can evolve independently and the effect of the bath is totally 
enjoyed by the friction coefficient rj. For the classical case, it is quite clear that at 
zero temperature the Brownian force can be neglected for the vanishing dissipation- 
fluctuation relation. According to the fluctuation-dissipation relation eq.(2.9), the 
classical statistical fluctuation of the bath is proportional to i]T . Only when T = 0, 
the system is isolated with an effective dissipation equation 



and the whole effect of bath on the system is only characterized through the friction 
constant 77. However, for the quantum case, the story is not so direct. Now, we 
consider the composite system C equals S plus B. Since C is closed its quantization 
process is well build in elementary quantum mechanics. 

Let us start with the exact solution to eqs.(2.2b) and (2.7) , 



q = E- rjq, 



(3.1) 



q(t)=Q(t)+X(t): 



Q(t)=a(t)q(0) + q(0)+g(t) 



N 



X(t)= £ *,■(*): 



(3.2) 



X j (t)=a k {t)x k (0)+/3 k (t)x k {0), 



9 



where 

1 - e~ vt E 

a(t)=t v = , g(t) = —(t -t v ) 

i] V 

C uj 
a k{t) = o , J 2 [ — sin[u k t\ (cos[r]t\ - 1) - rja(t)] 

C' i] oo k 
= L2,,a, [—{cos[u k t\ - 1) - -^-sin[rjt] + u k a(t)\. 

The canonical commutation relations at t = for the closed system of S plus B 

ih 

[g(0),g(0)] =ih, M0),i,(0)] = — (3.3) 

rrij 

quantize q(t) and q(t) at any instant t so that 

[q(t),p(t)]=ih, [x j (t),p j (t)]=ih (3.4) 

for the ordinary momentum- velocity relations 

p(t) = q(t),pj(t) = rrijXjit). 

Notice that the quantized coordinate q(t) is separated into two commuting part Q(t) 
and X(t) depending on the system and bath respectively. 

Now, we consider the quantum statistical problem of the dissipation and fluctua- 
tion for the above linear system. The quantum statistical average < >=< > qU antm 
over the bath is defined by 

<A>= i 13- 



for an observable A. Here, H B is the Hamiltonian for bath. A direct calculation gives 
the quantum fluctuation - dissipation relation 

< G(t) >= 0, 

D(t - Ot-xj = l< {G(t), G(t')} >= ^ r u jC oth{^p-)cos[uAt - t')]duj. (3.4) 
2 it Jo 2 
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At high temperature limit, ie., (3 — * or T — > oo , the above results approach the 
classical dissipation-fluctuation relation (2.9 ). In zero temperature limit, it becomes 



rqfo roc 

D(t — t') T _> = — / ujjCos[ujj(t — t')]dujj 
IX Jo 



- —Itm^ifi [ ^_ tl) - (M ^ )2 jj (3.5) 

This equation concludes that the Brownian force can hardly be neglected even at zero 
temperature for quantum case, but the corresponding fluctuation is proportional to 
rjh. So this is a quantum fluctuation. The effective coordinate Q(t) can be used to 
approach the physical coordinate q(t) only when this quantum fluctuation is neglected 
in certain sense. In next section we will clarify further the exact meaning of this 
argument in terms of the wavefunction description. 

When the quantum fluctuation can be ignored, the evolution of the dissipative 
system may be approached by the variable Q(t) independent of the bath. Then, the 
system is isolated from bath to produce a "closed " dynamics. Now, let us derive 
the effective Hamiltonian governing this kind of dynamics. Our derivation is valid 
for both the classical and quantum case 

By making an observation that the explicit expression (3.2) for Q(t) determines 
the commutator 

[Q(t),Q(t)]=ihe-*, (3.6) 
it is quite natural to define the canonical momentum-velocity relation 

P(t) = e r ' l Q(t) (3.7) 

for the effective variable Q so that the basic commutators for the Hamiltonian dy- 
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namics is 

[Q(t),P{t)]=ih. (3.8) 

Notice that the definition of the canonical momentum based on eq.(3.6) is not unique 
and the different definitions such as in ref.[18] may give different forms of the effective 
Hamiltonian. For the given Q(t) changing as eq.(3.2), one can explicitly obtain 
expressions of Q(t) and P(t) in terms of Q(t) and P(t) 

Q(t) = -Ee*Q(t), P(t) = e-*P(t), (3.9) 

Associated with Heisenberg equation 

Q(t) = ^-[Q(t),H e (t)], 

P(t) = ±:[P(t),H e (t)], 

the eqs.(3.9) lead to a simple system of partial differential equations about the effec- 
tive Hamiltonian H e (t) 

?M± = -Ee-»Q, d -Ml = e -*P(t). (3.10) 

A solution to eq.(3.10) determines an effective Hamiltonian H e (t) 

H e = H e (t) = ie^'P 2 - EQe rit + F(t). (3.11) 

up to an arbitrary function F(t) independent of P and Q. It can be regarded a 
generalization of CK Hamiltonian [14-16] of damped harmonic oscillator. 

In spirt of the above discussions, for an arbitrary potential V(Q), one guesses the 
general effective Hamiltonian for dissipation problem 

H E = H E (t) = l -e-*P 2 + V(Q)e r > t . (3.12) 
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Indeed, the Heisenberg equation of He = H E (t) with the basic commutator (3.8) can 
result in the dissipation equation 

Q(t) = - V Q(t) - ^) . (3.13) 

for arbitrary potential. 

In fact, the generalized CK Hamiltonian (3.12) can be directly derived in the 
explicit equations satisfied by Q(t) and P(t) as Heisenberg equation . It follows from 
eq.(3.13) that the commutator [Q(t),Q(t)] at time t satisfies an equation 

j t [Q(t),Q(t)] = -v[Q(t),Q(t)]. (3.14) 

Even in present of an arbitrary potential V(q), it still leads to the same commutator 

[Q(t),Q(t)]=ihe-*. (3.6') 

as that for case of harmonic oscillator and the case of constant external field. Then, 
eq.(3.6) suggests the same canonical momentum- velocity relation 

P(t) = e*Q(t) (3.7') 

as eq.(3.7) for the two above mentioned cases. Using the dissipative equation eq.(3.13) 
and eq.(3.7'), we have 

P( t) = -e-^M (3 . 15 ) 

Equations (3.15, 3.7') and the Heisenberg equation determine the equations of com- 
mutators about unknown-Hamiltonian H E 

[Q(t), H E (t)\ = ihe^P(t)/M, (3.16) 

[P(t),H E (t)] = -ihe-^?^-, 
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Obviously, the generalized CK Hamiltonian (3.12) is just a solution of above equa- 
tions. Therefore, we derive out the generalized CK Hamiltonian for the arbitrary 
potential V(q). Notice that the derivation of the generalized CK Hamiltonians (3.12) 
and its special case (3.11) here can also make sense for the classical case without 
Brownian motion so long as one use the commutators yinsteaded of the Poisson 
brackets. 

4. Meaning of Wave Function for Dissipation 

In the last section we derived the effective Hamiltonian in terms of the Heisenberg 
equations. It can serve as a starting point studying some dynamical problems of 
quantum dissipation under certain sense. To understand the problem completely, 
one must investigate the physical meaning of the (effective) wavefunction defined 
by the effective Hamiltonian (3.11) and thereby clarify in what sense this effective 
Hamiltonian can be used correctly. Thus, we turn to the Schroedinger picture and 
detail the main ideas and methods developed in ref. [1]. 

Let us first show that the wavefunction of the composite system C equals S plus 
B can be reduced to a direct product of wavefunctions in in two independent Hilber 
spaces. To this end, it is observed that the coordinate operator 

q(t) = Q(t)+X(t) 

of S is the sum of two commuting parts Q(t) and X(t). So the eigenfunction of q(t) 
with eigenvalue q is expressed as a direct product 

\q,t >= \Q,t > <8>|fi,i > <g>|f 2 ,f > ® \%N,t > (4.1) 
of the the eigenstates \Q, t > of Q(t) and t > of Xj(t) with the eigenvalues Q and 
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£j respectively. Notice that these eigenvalues satisfy 

N 

? = Q + Efc- ( 4 -2) 

3=0 

Due to eq.(4.2), for a given q, there exist many different sets of {£1, £2, 6v} cor- 
responding to q , that is to say, these eigenstates \q,t > are degenerate. So, the 
new notation \q, {£.,•},£ > with additional index {£,-} is needed, instead of \q,t >, to 
distinguish among the different degenerate eigenstate with the same eigenvalues q. 
Correspondingly, the following notation for arbitrary complex number £j 

\Q, {&} >= \Q > ®I6 > ®I6 > <E> ICjv > (4.3) 

is used to represent the degenerate eigenstate of q(0) with eigenvalue Q , which belong 
to the Hilbert space 

N 

V = V s <g> V B = V s <g> J] ® V i 
of the composite system C. Here, Vs and Vb = Ilf=i ®Vj is the Hillbert spaces of S 
and B respectively; Vj is the Hillbert space for the j'th oscillator in the bath B; \Q > 
is the eigenstates with eigenvalue Q and |£j > (e Vj) the eigenstates of Xj(0) with 
eigenvalue 

Let the composite system be initially in a product state 

N 

1^(0) >= 10 > >=i0>®n >> ( 4 - 4 ) 

at t=0 where \<p >, |VF > and \W,i > belong to V s , V B and Vj respectively. The cen- 
tral problem we should face is whether the product form similar to eq.(4.4), where 
the wavefunction is a direct product of two wavefunctions in two independent Hilbert 
spaces, is persevered in the Schrodinger evolution. The positive answer will implies 
that the system can be isolated from the bath by making use of an effective Hamil- 
tonian. To consider this problem, it is necessary to calculate the evolution operator 
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U(t) or its matrix elements. Because the coordinate operator q(t), for general value 
of t, can formally be regarded as an unitary transformation 

q(t) = U(ty q (0)U(t) 

of q(0), the eigenstate 

\q,{Q},t>=U(t)^\q,{Q}> (4.5) 

of q(t) with eigenvalue q can be constructed in terms of the evolution operator and the 
eigenstate \Q = q, {£,-} > of q(0) with the same eigenvalue q. Then, the coordinate 
component of the evolution state 

m) >=u(t)mo)> 

can be calculated as 

=< Qd&Mt) >=< QA^}\u(t)\m > 

= [< m\u(ty\Q,{Q >\* = [< m\Q,{Q,t >\* (4.6) 

according to the eigenstate \Q, > of q(t) for general t. In fact, this eigenstate 

can be directly solved in the coordinate representation with 



d ■ d 

9(0) = Q,q(0) = -ih—,Xj{Qi) = 0,^(0) = -ih/rrij — 



Obviously, the eigenstates \Q,t > of Q(t) and \£j,t > of Xj(t) with eigenvalues xij : 
MQ',t) =< Q'\Q,t >= e ^-W 2 HQ- 9 (t))Q'+HQ)] 



u^t) =< g|fc,f >= (4.7) 
are explicitly obtained from the differential equations 

HMO + Q' + a(t)]MQ', t) = QMQ', t) 
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(4.8) 



Now, we can write 



N N 



< q', {Q\q, {q, t >=< q'\q = q - t > n < t > 



3=1 3=1 



N 



=^E&^'*)n«&(4*), 



(4.9) 



which results in 



A? 



{s» = t< mh m, t >r = w{ q - e *r n o* (4.10) 



Notice that the variables £j are related to the bath, but they are not the coordinates 
Xj of bath. Therefore, applying the above analysis to a practical problem, one should 
distinguish between £j and Xj. 

It is observed from eq.(4.10) that , if the variable q is so highly excited that the 
Brownian motion contribution J2 £j is small enough in comparison with q in certain 
sense, the first factor in the right hand side of eq.(4.10) is approximately independent 
of £j. In this sense, the whole wavefunction for composite system C is factorized 
in to two independent parts belong to Vg and V B respectively. The first part W 
represents the wave function of the dissipative system isolated to evolve according 
to the effective Hamiltonian (3.11). Because of the Brownian motion, the physical 
variable q — Q + fluctuates about Q through the mean value of (J2Cj) 2 



where 



W(<7-E&) = JdQ'< m\Q' > <f> q -^(Q',t) 



N 



h 



<(Efc(*)) 2 >=E 



(\a 3 (t)\ 2 + ^\m)\ 2 )coth 



2KT 



(4.11) 



3=0 
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at temperature T. It is zero at t = and approach its final value 

<K« f =~)] 2 >=4^ +orcf O' (412) 

with low temperature limit in the time of the order of I/77. Obviously, whether the 
effective Hamiltonian can work well or not mainly depends on whether the values 
< (J2£,j(t)) 2 > can be neglected in practical problems or not. 

For the further study in the dynamics of the dissipative system, we need to cal- 
culate its propagator. There are usually two ways to do that in principle: 1). By 
solving the effective Schrodinger equation 

ih^Mt) >=H E (t)\9(t) >, 

directly for the obtained effective Hamiltonian H E . 2). By using reduced density 
matrix in terms of the path integral for the composite system C equals S plus B. 
Now, we will deal with this problem along a shorter road , as the third way in 
which we need not to know what is the effective Hamiltonian H E . Having known 
the time-evolution of the observable Q(t) explicitly, we can derive the propagator 
from Q(t) directly without use of the effective Hamiltonian H E . H E will thereby 
be derived in a purely quantum mechanical way. This approach not only avoids the 
complexity in calculations in the first two approaches, but also give us new insight 
for the understanding of quantum dissipation. 

According to the definition of propagator, we have 

G(q 2 , {£j,2}t 2 , qi, {£j,i}> *i) =< 02, {0,2}, t 2 \qi, {£,,1}, h > 

N 

= < Q2,t 2 \Ql,t 1 > J] <&,2,*2|&,1>*1>> 
j=l 
N 

= GiQ^Q^tJYlGji^hfaih). (4.13) 
18 



Because of the linearity of of the Heisenberg equation in the variables Q(t),Q(t) 
and £j(t),£j(t), the operators 0(t 2 )(0 = Q,Q,£j,£,j) at time t 2 must be a linear 
combination of the operators 0{t\) at ti. For example, 

Q(ti) = Q(h) + a(t u t 2 )P(t 2 ) + 6(fi, t 2 ) 

Q{h) = Q(h) + a(t 2 , + b(t 2 , h) (4.14) 

where 

Kt 2 , h) = g(t 2 ) - gfa) - Eaitje'^aih, h) 
a{t 2 ,t!) = -a(t u t 2 ) = a(t 2 ) - a(ti). 
Then, the definitions of eigenstates of Qi(t) 

Q(ti)\Qi,ti >= Qi\Qi,ti >,i = 1,2 

and its Qj, (j ^ i) -representation leads to the partial differential equations for the 
propagator 

d 

[Q 2 + a (t u t 2 ).ih— + b(t u t 2 )]G(Q 2 , h; QiM) = Q 2 G(Q 2 , t 2 ; Q u h) 

[Qi - a (t u t 2 ).ih-£- + b(t 2 , h)]G(Q 2 , t 2 ; Q u h)* = QiG(Q 2 , t 2 ; Q u h)* (4.15) 

Solving the above equations, we obtain the propagator 

1 



G(Q 2 , t 2 ; Qi, ti) 



xexp[ - - { l -Q\ + \Q\ - Q,Q 2 + b(t 2 , i^Qi + b(t u t 2 )Q 2 + 9(t)}] (4.16) 

where 6{t) is an arbitrary function of time independent of Qi, (i = 1, 2) Similarly, we 
can also calculate the factors Gj(£j,2, t 2 ; £j,i£i), but here we need not to explicitly write 
them out for our present purpose. For highly-excited q-system where the Brownian 
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motion is ignored, the first factor G(Q2, t 2 ; Qi, h) can be regarded as the effective 
propagators for the dissipative system. By taking into account that the propagator, 
in fact, is the certain matrix elements of the evolution operator Uq(t), i.e., 

< Q2\U Q (t)\Q 1 >=< Qi\U Q (tY\Q 2 >*=< Q 2 ,t\Q x >= G(Q 2 ,t;Q 1 ,0) 

the effective Hamiltonian (3.11) can also be derived again from 

H E = in?Mlu Q (t)-i (4.17) 

Notice that this derivation of the effective Hamiltonian here is purely quantum me- 
chanical. 

5. Spreading of Wave Packet Suppressed by Dissipation 

In this section the effective Hamiltonian (3.12) for the dissipative system is applied 
to study a quite simple dynamical problem: the motion of the wave packet of a 'free' 
(E=0) particle of mass M in one dimension. This study will show an interesting fact 
that the dissipation must suppress the spreading of the wave packet if the breadth 
of initial wave packet is so wide that the effect of Brownian motion can be ignored. 
Usually , without dissipation, a Gaussian wave packet will infinitely spreads into the 
whole space and the localization of wave will lose completely in evolution process. 
In future direction its breadth increases as t until infinity while its height decreases 
from its initial value to zero. Notice that the height and breadth of a wave packet are 
correlated through its normalization. However, for the present case with dissipation, 
there appears a quite different picture about the wave packet spreading. It will be 
proved that the final breadth and height have the finite limit value as t — > 00. In the 
following we use 1 and denote by x the operator Q approaching the physical 

coordinate q. 
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Starting with the effective Schrodinger equation 

governed by the effective Hamiltonian in absent of the external field, the evolution 
of wave function is generally expressed by 

q?( x , t)=J2< k\V(x, 0) > exp[ikx ^] (5.2) 



where 



2M 

is the energy of the momentum eigenstate \k >: 

< k\x >= -^=e ikx . 

Notice that, for the definition 

, x Mil- e-"*/ M ) 
f = a{t) = —± 

V 

we have t v — > t when 77/M — > 0. So one can regard as a ^-deformation of time 
t. Especially, t v approaches a limit — as t — > 00. This fact will enjoys the physical 
features of wave packet spreading in present of dissipation. 
Take the Gaussian wave packet 

0) =< x|*(0) >= L^_ e **o*-£ ( 5 3 ) 

[27T(i 2 ] 4 

as an initial state, in which 

< *(0)|x|*(0 >= 0, < tf(0)|P|tf (0 >= hko, 



< Ax >= J< ty(0)\x 2 \^(0) > - < *(0)|a;|*(0) > 2 = d. (5.4) 
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The above equation (5.4) shows that the wave packet is centered at x=0 and has an 
average momentum hko. Its breadth is d. According to the wave equation (5.1), the 
Fourier transformation 



1 

/f?r 
2_f 

7T 

for the initial wave packet (5.3) determines the wavefunction at t 



j roo 

*(x,0) = -= i){k)e ikx dk 
\/2ir J co 

^(jfe) = _-l e - d2 ( fc - fc o) 2 



( ' )_ (2ny/^d + U v h/(2Md) eXP[ ^ U ^ M) d^(t v h/2Mdf ] ^ 
To understand the physical meaning represented by the wavefunction (5.5), we 
write down the corresponding position probability: 

\*(x t)\ 2 = 1 exv\ ( x -k Q t v h/Mf 

1 1 ' )l ^TM^S 2[^ + M ) V (2M^] J - ^ 

The above formula tells us that the motion of the wave packet in dissipation process 

starts with initial velocity 

Vo =< p/M >= hk /M 

and its center is initially at the position x=0. Subjected to time evolution, the center 
of the Gaussian wave packet will stop at a limit position 

hko 

% limit • 

V 

as t — > oo. In this process, the velocity of the center 

A fit. 

v(t) = j t < y(x,t)\x\^(x,t) >= |e-"/ A ' 

decrease from ^ to zero. The above fact means that the motion of the center of 
wave packet is as the same as the that of a dissipative classical particle. However, a 
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purely quantum picture is manifested by the finite change of its breadth 

B{t) = yjd? + (kt v /2Md) 2 
from d at t=0 to a limit value 

B limit = y/d? + (h/2 V d) 2 

as t — > oo. Both the position at which the center of wave packet finally stop and 
its final breadth are independent of the mass M! It defines the limit shape of wave 
packet at which the spreading wave packet finally takes. These physical features are 
illustrated by the figures 1 and 2. 

It is finally pointed out that this suppressing of wave packet spreading by dissipa- 
tion possibly provides a mechanism to localize a quantum particle. 
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Figure Captions 
Figure 1 

The Motion of Wave Packet in Present of Dissipation : F(x, t v ) = 
\^(x,t)\ 2 . The dissipation not only restrict the motion of the center of the 
Gaussian wave packet like a classical particle, but also suppress its spread- 
ing so that it takes a limit Gaussian wave packet with finite breadth and 
height as t — > oo. 

Figure 2 

The projection of Figure 1 on the |^| 2 — x Plane. It actually represent 
the figures F(x,t v ) = \^{x,t)\ 2 at t=0, t=T, t=2T, t=3T ,...,t=NT. When 
iV — > oo, there is a lower Gaussian wave packet \^(x : t = oo)| 2 = F(x, 1/rj). 
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